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Abstract 

Motivated by Polychronakos’ discovery that solitons exist in the hydrodynamic equations of 
continuum version of the Calogero model, we seek solitons in the classical dynamics of a continuum 
version of the Haldane-Shastry spin chain. We have obtained analytic multi-lump solitary wave 
solutions for our spin-field equation, and these solutions possess interesting topological features. We 
have performed numerical collision experiments showing that these solitary waves survive collisions, 
and thus suggest the existence of true multi-soliton solutions. 
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I. INTRODUCTION AND MOTIVATION 


Solitons are localized wave packets that survive unchanged through collisions. Since their 
accidental discovery as water waves described by the KdV equation (See [T] for an account 
of this), soliton solutions have been found in many systems of partial differential equations. 
They have also been observed in many physical systems — for example, nonlinear optics[2l|3], 
matter wave solitons in Bose-Einstein condensation|3], and vortex rings in ferromagnetic 
materials [5H7]. 

The soliton property is closely related to the integrability of the underlying system of dif¬ 
ferential equations. Classical integrability ensures that there are as many Poisson-commuting 
integrals of motion as degrees of freedom, and the associated conservation laws suppress the 
available post-collision phase space to the extent that the only possible outgoing solutions 
are non-diffractive and just a permutation of the incoming solitary waves. 

It is natural to describe a quantum integrable system as one with as many mutually com¬ 
muting local operators as degrees of freedom. This is not necessarily the most useful dehni- 
tion as knowing such operators does not always help solve the system. Sutherland therefore 
dehnes a quantum integrable system as one that supports non-diffractive scattering [8]. A 
beautiful example is provided by the Calogero-Sutherland models which are both classi¬ 
cal and quantum integrable [9]. It is possible to take a continuum limit of the classical 
Calogero model and the resulting hydrodynamic equations of motion possess multi-soliton 
solutionsfini E]. This remarkable result suggests that the continuum model remains classi¬ 
cally integrable. 

With Calogero hydrodynamics as motivation, we here construct a classical version of 
Haldane-Shastry (HS) spin chain [T2IITT] . The HS spin chain is the inhnite mass limit of the 
quantum integrable spin-Calogero model [TSHTT] and is also quantum integrable [18]. 

To obtain a classical continuum version of the spin chain, we hrst interpret the spin-^ ex¬ 
change term as ferromagnetic interaction between the spins, and then replace spin-(j = 1/2) 
by a large enough value of j so that the dynamics becomes classical. The degrees of freedom 
of our classical model are therefore unit vectors mj at equally spaced lattice sites. We take 
the interactions to be ferromagnetic so that we can anticipate a continuum limit in which 
rrij —)■ m(x) with m(x) being a smooth function. The resulting equation of motion for m(a;) 
is then a non-local generalization of the (known to be integrable) Landau-Lifshitz equation 
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where the second derivative with respect to x is replaced by the derivative of a Hilbert trans¬ 
form in a manner reminiscent of Benjamin-Ono equation [121120] (itself integrable equation 
having known multi-soliton solutions [21] )• We will see that our generalized Landau-Lifshitz 
equation possesses both interesting analytic solutions and numerically-obtained multi-soliton 
solutions. 

The paper is organized as follows. In section [TT| we construct our continuum model 
and derive the classical equation of motion. In section mil we introduce numerical and 


analytical methods, and use them to solve the single-speed sector. In section [rV] we presents 
numerical collision experiments of multiple solitons. Section |V| is the discussion of the 


solutions, especially their topological features. [V^ concludes our results. Some technical 
details are presented in the appendices. 


II. THE CLASSICAL AND CONTINUUM HAMILTONIAN 


A. The Continuum Hamiltonian 


The ferromagnetic version of the original Haldane-Shastry model has Hamiltonian 


^HS = 


E 


1 (T^ * (Tj 

{xi - XjY ■ 


( 1 ) 


Here cr^ = Cy, az)i are Pauli matrices that act on the Hilbert space of a spin at position 
Xi and we have taken 1 — cr^ ■ (jj in the interaction so that the energy would be zero if all 
spins were parallel. 

The spin-1/2 dynamics is very quantum mechanical. There are two routes to modifying 
the system so as to obtain a classical limit. One is to replace the SU(2) spin group with 
SU(iV), while preserving the interpretation as a spin exchange interaction. Thus 


Pij - 2 ^^ + 


( 2 ) 


where A = (Ai,... Xn^-i) are generators of SU(A) normalized so that tr(AaAfe) = 6ab- The 
classical approximation then has the A* taking values in a suitable co-adjoint orbit. 

The other, which is the route we will take, is to regard the (Ti -cTj part of the exchange term 
as a ferromagnetic interaction between spin-(j = 1/2) particles and then take j large enough 
for the spins to become classical. The resulting interaction term can no longer be interpreted 
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as a spin exchange, as the general spin-j exchange term is a higher-order polynomial in the 
spin operators. 

Our original HS Hamiltonian has therefore been replaced by 


^HS ^ H, 


classical — 


1 — rrij ■ nij 

^ {Xi - XjY 

1 <J ^ ' 


(3) 


We place the spin chain on a Id lattice with spacing a, so a ^ is the density p and we can 
write 


H, 


' m, — m. 


classical 


2 ^ 

P Z. Z. 7737 


J = -OQ 


jY 


(4) 


We now wish to replace the nij by the smooth function m(x) and the sums over i and j by 
integrals. 

There are two ways to approximate sums by integrals, one is to use real-space Euler- 
Maclaurin summation formula, and the second makes use of Fourier interpolation in mo¬ 
mentum space. The two methods give the same result (see appendix]^: for slowly varying 


mi 


H, 


classical 


Ho = p^ 


dx{'Km- 


Here m-^(x) is the Hilbert transform of m(a:) dehned by the principal-part integral 




m(0 


di. 


On using the identity 


d 


r izM, 

dx \n X - 


TT x-i 

^ mi{x) - mi{0 


di\=- 


TT 


{x - 0^ 




(5) 


( 6 ) 


(7) 


together with |mp = 1, we see that the hrst term in the integral is the naive continuum limit 
where we simply replace each sum by an integration. The second term is a correction that 
we have kept because such corrections play a vital role in preserving the soliton property in 
the classical hydrodynamics of the Calogero Sutherland models. It is informative to look at 
the role of this correction in momentum space. We have 


Ho = 


= P 


' —OO 
COO 


^ \ —isgn.{k) ■ ik ■ n — ^ |m(A:)| 


dk 


27r 


— < 7r\k\ - k^ > |m(/c)f. 


( 8 ) 


We see that the correction term leads to our continuum approximation to the manifestly- 
positive discrete-model energy density being in danger of becoming negative when it has large 
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k component. Of course the continuum limit is supposed to be smooth and so large k are not 
supposed to appear. In particular, values of |fc| > tt are meaningless since they alias to k — 
27r. Nonetheless, using this Hamiltonian leads to instabilities in the numerical simulations. 


The evolution and gradient-descent mentioned in section all develop unwanted anti¬ 
ferromagnetic {k = vr) oscillations. We therefore tentatively discard the double derivative 
term — in the Hamiltonian. 

Apart from the numerical instability, there is a legitimate reason for the discard. In 
momentum space, the Hilbert transformed term scales as |fc| and the double derivative term 
scales like /c^. We only care about the physics in large distance, so k is very close to zero. 
The Hilbert transformed term is thus far larger than the double derivative term and, in the 
hydrodynamic limit it is legitimate to neglect the later. The Hamiltonian we actually use is 
therefore 

1 

dx m ■ (9) 
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The constants and tt have also been dropped as they only affect the time scale. 


B. The Equation of Motion 


The classical motion of a single unit-vector spin m = (mi,m 2 ,m 3 ) is derived from the 
Poisson bracket 

= eijkmk. ( 10 ) 

We extend this bracket to functionals ^[m], G[m] of a continuous spin-held m(a;) by setting 

{mi{x),mj{x')} = eijkmk S{x - x'), (11) 


and hence 


{F, G} — eijk 


6F SG 


/_oo Smi{x) 6mj{x) 
The time evolution of m for a Hamiltonian id [ml is then 


mk{x)dx. 


(9m 




6H 


= -— X m. 

om 


For example 


1 

HLL[m] = - dx 

^ J —00 


( 12 ) 


(13) 


(14) 
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gives Landau-Lifshitz equation (LLE) [23l 


9m 9^m 

= m X 


(15) 


dt *** dx"^ ’ 

whose soliton solutions are connected to helical curve motion |2SH2H] and have been exten¬ 
sively studied [2^ 1221 - 152] . Many of the methods used for the LLE turn out to be applicable 
to our problem. 

For our classical continuum version of the Haldane-Shastry spin chain, the equation of 
motion becomes 


9 m 


= 9 a;m^ X m 


( 16 ) 


Because differentiation with respect to x commutes with taking a Hilbert transform, equation 

( 17 ) 


(21) can also be written as 


9 m f 9 m \ 

^ J. ^ 


H 


C. Conserved Quantities 


Integrable systems with continuous degrees of freedom process inhnitely many conserved 
quantities that in some cases enable us to construct multi-soliton solutions [33]. Here we 
will only discuss the more obvious conserved quantities associated with the systems’ global 
symmetries. The resulting constants of the motion will not only be useful in the numerical 
calculation but also characterize the general physical properties of our model. 

The Hamiltonian is time-translation invariant, hence the Hamiltonian itself as the time- 
translation generator is conserved. Also the invariance of Hamiltonian under global rotation 
ensures that each component i = 1, 2, 3 of the total spin 

/ OD 

mi{x) dx 

-OO 

is conserved. 

A space-translation invariant continuous spin chain should possess a conserved momen¬ 
tum, but a rotationally invariant expression for the generator of space translations was long 
missing. This issue was elucidated by Haldane who showed [3l| that the quantum-mechanical 
generator of hnite translations is the exponential of a Wess-Zumino term and is only well 
defined for discrete translations past an integer number of spin-l/2’s. If a local expression 
is required it necessarily involves a monopole gauge held in spin-space and the position of 
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the Dirac string breaks rotational invariance. A similar issue arises in the purely classical 
chain. It is nonetheless useful in numerical calculations to introduce the quantity 


P= Aim) ■ c?T.m dx 


(18) 


where A is a monopole gauge held. Provided m(a;) avoids the Dirac string at the south pole 
m = —z, we can take 

. (— z) X m , , 

A = A-A-r. 19 

If we treat m(a;) as a curve on the unit sphere parameterized by x, then P is the area of the 
region enclosed by that curve on the side without the south pole. 

We then hnd that dxiti = {m, P} for this restricted class of conhgurations. Because of 
the 47r ambiguity in dehning the “area enclosed” by P it is only the quantities 


U = e”'’, a/2 6 Z, 


( 20 ) 


that are well-dehned for general conhgurations. 


III. SINGLE-SPEED SOLUTIONS 


We begin by seeking a single-soliton solution to (16) that moves at a hxed velocity v, and 


so is of the form m(a;, t) = m(a; — vt). The equation of motion then reduces to a non-linear 
ordinary diherential equation 

rJim 

( 21 ) 


dm dm-^ 

V— — = m X —-— 
dx dx 


or, equivalently, 


dm f dm\ 


( 22 ) 


n 


It is not obvious that equation (21) possesses any interesting soliton-like solutions, so we 
hrst made a numerical search. 

We followed the strategy used by Tjon and Wright [22] to hnd the Landau-Lifshitz soliton. 


We note that equation (21) can be rewritten in term of the Poisson bracket as 


{dt + vdx)m = {m, H + vP} = 0, 


(23) 


and also note this same Poisson bracket equation gives the stationary points of the functional 
d/[m] +uP[m] subject to the variation on sphere. This because a variation a functional /[m] 
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under 5m = m x r; is given by 
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j {m, 1} ■ rj dx. 


(24) 


Thus requiring 5/ = 0 for all r) is the same as requiring {m,/} = 0. Consequently all the 
single-speed solutions are stationary points of /o[ni] = if[m] -|- nP[m]. 

Numerically accessible stationary points should be extrema[3S], so it is convenient to 
modify the functional Jq and include a penalty term that makes the functional positive 
dehnite. If we take /[m] = H + ci(P — Pq)^, at least one stationary conhguration can be 
found by minimizing the functional and so solving 


{m, 1} = {dt + 2ci(P - Po)9^)m = 0. 


(25) 


As momentum is a conserved quantity, equation (25) is identical to equation (21) once the 
velocity is identihed as 

n = 2ci(P-Po). (26) 


A naive gradient descent to a minimum would take, 


nit+i = mt + heffAt, 


leff — 


5/(m) 

5m 


(27) 


where the gradient heg is the effective field(mean field) produced by all the other spins and 
any external held. However, this strategy does not preserve the unit length of the spin held. 

Instead, we project the gradient in the direction perpendicular to m in each updating 
step[36l EZ] and take 

mt+At = [heflf - (heff ■ m)m] At. (28) 

In the At —>■ 0 limit, this becomes the Gilbert damping equation 

ctm 


dt 


= — m X (m X h 


effj 


(29) 


Both the Gilbert damping equation (29) and our generalized Landau-Lifshitz equation (16) 


are time evolution equations of the type, 

ttm 


dt 


= ca m X m 


( 30 ) 


We wish to solve them accurately together while preserving the condition |mp = 1. 









Among many numerical schemes proposed for Landau-Lifshitz dynamics [371 H01 - H2] we 


found a mid-point finite difference to be the most suitable. This sets 


nif+Af - mt 

At 


^K+iAt] X 


mt+At + nit 
2 


(31) 


and automatically preserves the length. In spite of being an implicit method, it can be made 
explicit by use of a predictor-corrector procedure. 

We found that the relaxation process works cleanly and nearly all smooth initial condi¬ 
tions converged to a single-speed solution. The strength of the residual effective held |heff| 
is a measure of the hnal numerical error and is of order 10“®. We can also time-evolve the 


numerical hnal solution under equation (16) and conhrm that it moves at the calculated 


constant velocity with an error of order 10 
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FIG. 1: A typical set of numerical results obtained by variational method. a),b) and c) are the three 
components of the single soliton solutions respectively, d) shows the trajectory of the minimal of ttIz in 
X G [—0.5, 0.5]. A least square fit of line gives ^measure = 0.5177, while inserting the numerical value of P 


and Po for this configuration into equation (26) gives v = 0.5176. 
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Figure [^demonstrates a typical set of data and fitting result. The measured velocity in the 


evolution agrees with the coefficient in equation (26), and so conhrms the existence of the 
single soliton solution. 

Having reassured ourselves that there was at least one stable single soliton solution, we 


sought an analytic solution to (21) by exploiting the fact that if a smooth real function u{x) 


obeys some mild condition — such as that it lies in some > 1 for example — and 

set v{x) = U'n, then the function. 




u(x) + iv(x] 


dx 


X — z 


(32) 


is analytic in the upper half plane H, tends to zero at inhnity there, and has boundary value 
on the real axis 

f{x) = u{x) + iv{x). (33) 

We therefore introduce a complex vector field M(2;) whose components are analytic functions 
on the upper half plane and whose boundary values on the real axis are 


M(a:) = m(a:) + im-^(x). 


(34) 


From the traveling-wave equation (21), the real-axis normalization |my = 1, and various 


Hilbert transform identities, we can deduce several results 


1. A decomposition of dx'm.'H'- 


dxva.'n = —um x dxva. + \/Y--^\dxVa\m.. x G 


(36) 


2. Two orthogonality relations: 


■ dxT^n.q-i = 0, d^m. ■ = 0. a: G 


(36) 


3. An analytic-function version of the orthogonality relations: 


dM ■ dM = 0, ■ d^M = 0 X G e 


(37) 


For proofs of these results see appendix [^ 


The hrst identity (35) tells us that a solution can only exist when the velocity is less than 


unity. Further, the complex orthogonal relation is very powerful. It simplihes the problem 
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by trading the non-local Hilbert transform for the analytic function M. Although initially 


derived by assuming that 2 ; lies on the real axis, both left-hand sides of (37) are analytic 
functions of z, and so both equations must also hold in the entire upper half plane H. We 


can now solve the extended (37) by expressing M 3 in terms of Mi and M 2 (Surprisingly, the 
solution implies = 0 for all n > 2 ), 


- {8^2?? = 0 

In -Tr—r—j = 0 

dzM2 

dzM 2 = C 2 dzMi similarly dzM^ = c^dzMi 


1 -|- C2 "h C3 — 0. 


(38) 

(39) 

(40) 

(41) 


The result is that the three components are proportional to each other. 

After using the global rotational symmetry, the general solution M(z) can be parameter¬ 
ized by a single analytic function g[z) and velocity v = cos^, as 

M(z) = l),isin0(^- l),l + sin^0(^- 1)^ (42) 

The point-wise constraint |m(a;)p = 1 that must hold on the real axis, together with bound¬ 
ary conditions (values of m(±cx))) further dictate that 


g{x)g[x) = 1 g{x = ±cx)) = 1. 


(43) 


A general solution can now be obtained by applying two-dimensional potential the¬ 
ory. The analyticity of the real and imaginary parts of the meromorphic functions 
<h(z) = —^\Yig[z) = u + iv implies that they satisfy the Poisson equation. They can 
therefore be regarded as the electrostatic potentials produced by charges located at zeros of 
g{z). The condition g[x)g{x) = 1 sets the boundary conditions uiy = 0) = 0. Hence u{x, y) 
is just the potential of a semi-inhnite metallic slab and obeys 


-VM 


H/{zeros of g} 


= 0 


u{y = 0 ) = 0 . 


(44) 


By the method of images, we can remove the boundary condition at the price of placing 
an image charge at Zi in the lower half plane for each charge located at 7^. A general 
hnite-charge solution is therefore 


N 

9{z)= n 

i=l 


z - Zj 
z- Zi 


( 45 ) 
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For each such meromorphic g{z), a spin held solution moving with velocity v = cos 6 is then 


nisoi = (-sin(26')(5(/j - l)/2, - sin 1 + sin^ 9{gR - 1)). (46) 


It is remarkable that these single-speed solutions possess such a rich structure! 

The reasoning leading to this general single-speed solution is quite subtle, so as a reality 
check we have conhrmed — both by substituting the proposed solution into the equation 
of motion and numerically — that we have indeed found multi-lump single-speed solutions 


to (21). Because of the periodic boundary conditions used in our numerical work, only g{z) 
with periodic poles are candidates for comparison. The simplest case is the periodic version 
of single pole solution. 


9{z) = JJ 


z — n — ia sin 7i{z — ia) 
z — n + ia sin 7i(z + ia) 


(47) 


which reduces to when z a. Setting a and velocity v as htting parameters, the 
least square ht shows that the numerical and analytical results match perfectly with only a 
residue of order 10“®, which is roughly the order of numerical error. The htting velocity for 
the conhguration shown in £gure[T]is 0.5176, the difference is less than 10“®. 


IV. NUMERICAL MULTI-SOLITON SOLUTIONS 


The existence of single-speed multi-lump solutions suggests that there will be multi-soliton 
solutions that move at different speeds. We have not, however, been able to generalize the 


pole ansatz (45) so as to decouple the speeds and so find such solutions analytically. 


We therefore arranged several single-speed single-lump solutions some distance apart in 
the hope of approximating the initial conditions of a multi-soliton collision. We then evolved 
the spin held according to (16). Typical soliton interaction behaviors are shown in hgurej^ 
and 1^ The two solitons not only retain their initial prohles, but also appear to have zero 
time lag in the asymptotic region, as if the collision never occurred. 


Collisions of three solitons are shown in hgure and Their behaviors are qualitatively 
the same as the superposition of three 2-body collisions. 
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FIG. 2: Collisions of two single solitons. The 
displayed amplitude is 1 — for better 
visualization. The velocities for the right-going and 
left-going solitons are 0.9 and —0.8 respectively. 


FIG. 3: Color map representation of the left figure. 
The larger amplitude soliton suffers less than the 
smaller one, so the time lag and advance are 
insignificant except for the region close to the 
center. The smaller starts to deviate earlier, but in 
the end either have net time lags. 


V. DISCUSSION 


Through numerical and analytical calculations, we obtained a large family of single-speed 


multi-lump solutions to the generalized Landau-Lifshitz equation (16). The speed is always 


less than 1 and serves as a parameter that controls the amplitude. The shape of each lump 
(dehned as the full width at half maximum for each Lorentzian solution ) can be different 
depending upon the imaginary part of the pole. 

One curious feature of the single-speed solutions is that the total energy and momentum 
are completely ignorant about the positions of the poles in g{z). 

Consider the energy for a general spin held where M = (/i, / 2 ,1 -|- /a) with three possibly 
different analytic functions as its components. Each component fi maps the upper half-plane 
El to a region Di. One can show that, 

^ ^ . r rt Ir/f- 

(48) 


E = i / M ■ d.Mdz = 


2i 


2i 


^=l JdD, 


hdh = Y. 

i=l 


f dfjdfi 

Id, 


Hence the energy is the sum of the oriented areas swept out by the three analytic functions. 
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FIG. 4: Collisions of three single solitons. The FIG. 5: Color map representation of the left figure, 
velocities are 0.8, —0.25 and —0.85 from left to right 
in the initial state. The three-soliton collision in 
this figure can be viewed as a ’’superposition” of 
several 2-body collisions. 

For the single-speed solutions this gives 

E = 28111^ e f = 2TTNsm^e = 2ttN{1-v^), (49) 

J unit disk 

which is proportional to the winding number N of g. 

Similarly the momentum P is the solid angle enclosed by m(x)[3l]. Hence it should 
depend only on the global property of the spin held. For the single-speed solution, it is also 
proportional to N, 

P = 27rAf(l — cos 6*) = 2tt N{y — 1) (50) 

The mysteries of these neat results can partly be elucidated by the geometry of the 
parameterization ( [42| . By taking a scalar product, we hud the projection of m on the titled 
unit vector n = (sin 6*, 0, cos 6*) is the constant velocity cosO. So the trajectory of the tip of 
m on a unit sphere is in fact a small circle. Thus a more nature parameterization is to take 
n as the z axis, and m to be 

m = (gfjjsin 6^, sin^, cos 0 j (51) 
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where g = gR + igi is the same analytic function dehned in equation (45). The two param- 
eterizations of m are related through a rotation about the y axis and thus are equivalent. 
Nonetheless the geometric meaning is clearer in equation (51). As shown in figure]^ for a 



FIG. 6: The single-speed solution m(a;) viewed 
as a parametric curve of x on unit sphere. When 
X runs along the real axis, the tip of m traces out 
a circle on constant latitude z = v = cos 9 and 
repeats N times for an A^-lump solution. 


Af-lump single-speed solution, the tip of m traces out a small circle on the constant latitude 
plane z = v = cos 6 and repeats with possibly different paces N times on this circle. The 


momentum 2Af7r(l — cos 6) calculated in equation (50) is actually the area of N spherical 
caps. Furthermore, the image of the analytic function Mj, on the complex plane is a circle 
congruent to small circle in figure Consequently, the area swiped by is area of N of 
those small circles, which is Nnsin^O. The same is true for My. Taking into account that 
Mz is constant, the total oriented area is thus 2tiN sin^ 6 as we calculated before in the total 


energy equation (49). 

It now becomes apparent that all the dynamics of single-speed solution is constrained on 
a circle, which collapses the possibly existing multi-soliton space to the boring single-speed 
solution space. Generalizing the U{1) group element g{on the real axis) in the single¬ 
speed solutions to the one in a larger group(e.g. SU{2)) is necessary to accommodate the 
numerically-found multi-solitons. 
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VI. SUMMARY AND CONCLUSIONS 


We have proposed a classical version of a continuum Haldane-Shastry spin model with 
the Hamiltonian H = dxm ■ The non-local Hilbert term is a consequence the of 

long-range interactions in the original quantum model, and the equation of motion contains 
a Hilbert transform that is similar to that in the integrable Benjamin-Ono equation. Moti¬ 
vated by the results of Polychronakos on the continuum Calogero-Sutherland model and of 
Abanov and Wiegmann|43j on the Benjamin-Ono equation, we conjectured that the model 
is integrable and support multi-soliton solutions. 

We numerically sought for and found single-soliton solutions that move with constant 
velocity |n| < 1. We then found an analytic form for a large class of single-speed solutions 
by rewriting the equations in terms of analytic functions on the upper half-plane H. These 
solutions contain multiple lumps with possibly different widths. We found that the energy 
and momentum of a V-lump solution are proportional to the topological invariant winding 
number N. From the geometric point of view, it is the consequence of the circular trajectory 
shown in hgure 

We also performed numerical soliton-soliton collision experiments and found that the 
single-lump solitons survive multiple collisions unscathed. There appears to be no asymp¬ 
totic time lag after the collisions, and this, and the role of the winding number N, suggests 
that there should be some clever transformation to a system of non-interacting particles. At 
the moment we have no idea of how to hnd this transformation. 
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Appendix A: Two methods to calculate hydrodynamic limit 


The double summation we are going to do is 

i/ = p2 ^ 




■jy 


^ nii ■ Si 


(Al) 
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where Si = } is the inner sum. The total sum and the inner sum Sj are both 

j=—oo *•* 

absolutely convergent, so it is legitimate to reorder and do the inner sum first. 

Our first approach is to use a special Euler-Maclaurin summation technique introduced 
in reference na. Euler-Maclaurin formula converts a discrete sum into an integral of the 
interpolation function f{x) 


n=a+l 


fin)= I fix)dx + 


k=0 


k\ 


-\- Rpi Qj^b ^ Z ,CL <c. b. 


(A2) 


If the interpolation function is smooth and has exponential type less than 27r, then the only 
corrections are derivative at boundaryies and the remainder term 


Rp = (- 1 )^ / -jBp{x - [x\)f^P\x)dx, 
J a P- 


(A3) 


where Bernoulli numbers and Bernoulli polynomials Bp{x) are accompany with the pth 
derivative. 

Due to the singularity at z = j, one can not naively apply it to the inner sum Sj . However 
it is possible to introduce a counter term as in |13] to remove the singularity without affecting 
the result. The regulated summand is fjj = ~ ^ so that 

j=N^ 

Si = lim } fjj = lim 

TV-s-oo N^oc 

j=-N 


-V 


dz/f,;,. 


>-N 


(A4) 




limfj^-F lim 


k=l 


-N 


+ Rn- 


(AS) 


We only focus on the soliton solutions, which vanishes asymptotically f^j\±oo) = 0. Hence 
the derivative corrections at ±N go to zero for hxed order p. The remainder term has a 
bound 

(A6) 

for the finite domain summation, the remainder term can be arbitrarily small by taking a 
large enough order p. For the inhnite domain, in most cases the remainder term can also be 
neglected for functions whose exponential type is less than 27r( It is however not necessarily 
zero, see an example of applying this formula to the Id monoatomic gas partition function 
in chapter 7 of the book|45]). 
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Granted that we can neglect the remainder term, the discrete sum has only two terms 
left 

Si = (A7) 

where the Hilbert transform comes from the integral, the double derivative comes from fj*. 
Completing the sum over i gives the total energy. 


H = 


^ poo I' 

Si ■ nii = y d/ij -m(/i)m"(/i) 


fOO 

= P^ d/i[—-m'^ + 7rm ■ m^]. 
J — OO ^ 


(AlO) 


An alternative way to do the sum is the Fourier interpolation. The advantage of this par¬ 
ticular interpolation is that the discrete sum is automatically equal to the integral, because 
the zeroth component of Fourier series is dehned to be the sum! 

The interpolation starts with expanding a periodic function f{n) by its discrete Fourier 


/(”) = fkt 


(All) 


2ir 2 ^ 


and then extending the dehnition of Fourier series for integer number n to real number x 


f(^) = Y 


(A12) 


What immediately follows is that 


N 1 

= / 

n=l >^0 


f{x)dx. 


(A13) 


Now take any component of m to be a periodic function f{n), the corresponding compo¬ 
nent in S; is then 


y- 9(1)-gjj) 


E /-“E^V 


jt—f+i j—“ 


E A^^'E 


' 1 — e* 


(A14) 


(A15) 


1 ^—n+i 

27r 2 
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Identity (A.4) in the reference |16] provides a closed form expression for the summation over 
n. In the notations of this appendix, the identity is 


Li„(e*'=) + (-l)"Li„(e-*'^) = 


n\ 


'2ti' 


where Li„(a;) is the polylogarithm function 

which is equivalent to 


n=l 


E' 


Akn 






We need the case of s = 2, so that for 0 < \k\ <2 ti 


_. Akn 




Therefore 


E' 


j=-oo 


9{i) - gjj) 
ij -1? 


w 

2 




27r 2 ^ 


and the inner sum becomes 


OO / 7 \ / * \ 

9{l) - dU) 1 ^// , 

E ^j_i)2 =of +^^(/)L=r 


j=-oo 


(A16) 


(A17) 


(A18) 


(A19) 


(A20) 


(A21) 


The result of Si agrees with (A7). Moreover, there is no approximation involved in 


replacing the outer sum by integral for Fourier interpolated functions. Therefore the total 


energy expression H in equation (AlO) becomes exact in this method. 


Appendix B: Equivalent forms of traveling wave equation 


This appendix presents the algebras to deduce several equivalent forms of the traveling 


wave equation (21). 


We begin by showing equation (21) is equivalent to a decomposition 


= —nm X dx^i^ + \/w4F|c?3,m|m. 


(Bl) 
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Assuming is nowhere zero, the set {m, m x c?a;m} forms an orthogonal basis. 
So dxm.H can be expanded as a linear superposition of them. Taking the scalar product of 


the traveling wave equation (21) with we hnd = 0. Hence there is no 


c? 2 ;m component in the decomposition of dx^-'u 


dxm.'H = «m X ^xy^^ + /dm. 


(B2) 


Taking a cross product of rtB2h with m, the consistency with equation (21) requires a = —v. 


To determine /d, we expand ^{dxm ■ ^^m) by the Hilbert transform identity for the product 
'Hifg) = fn9 + fgn + Wngn)^ to get 


nidxni ■ dxui - OxUiH ■ dxmn) = 0 . 


(B3) 


Whatever in the braces is a constant, which is set to zero by boundary conditions. So the 
two vectors have equal lengths 


dxin ■ dxin = dxinn ■ dxinn, 


(B4) 


The “equal-length” condition together by squaring equation (B2) allows us to compute jS 
|c?a,m^| = -I-/d^ 13 = ±\/l — v‘^\dx'm\. 


(B5) 


The positive dehniteness of the total energy m ■ dxin-fi = f3dx > 0 picks out the 


positive (3 branch. We arrive at the decomposition (Bl). Since the converse is trivial, this 


completes the proof of their equivalence. A byproduct is the that the speed is always less 
than 1. 

Furthermore, the decomposition we just derived is also equivalent to the two orthogonality 
relations 

dxin ■ dxin-H = 0 = 0. (B6) 


We have proved the hrst in equation (B4). The second can be verihed by differentiating the 


decomposition relation (Bl) and doing a scalar product 


dim ■ d'lm'n = [ — vm x d^m + \fl^^^dx{\dxm\m)\ ■ d\ 


m 


= vT 


(c?a,m ■ dim) 


\dxin 

[(m ■ dim) |9a;mp] 


= 0 . 


(B7) 

(B8) 

(B9) 

(BIO) 
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So the second derivative terms are also perpendicular to each other. Similar to the derivation 
in equation (B3), we can even prove an analogous “equal-length” relation 


dim ■ dim = dlm^ ■ dlm^. 


(Bll) 


On the other hand, given these two orthogonality relations, we can also deduce equa¬ 


tion (21). The hrst orthogonality relation tells us that d^m^^ can only have the following 


decomposition (note that a may be position dependent) 

d^rny^ = am x dxm ± a /1 — a^\dxm\m. 


(B12) 


From our calculation above, ■ dim = 0 if a is constant. So the scalar product must 

be proportional to the derivative of a 


dim ■ dlm-u = dxCi{m x dxm =)= 


a 


\/l — 


a^ 


;|c?a;m|m) ■ dim 


dxO. 


a 


dim ■ dlmn =F 


dxCx 


a 


\/r 


a^ 


;|9a:ni|m ■ 


= ±- 


dxd 


a\/T— 

the second orthogonality relation restricts a to be a constant. 

In summary, these derivations rewrite the traveling wave equation in completely differ¬ 
ent forms without losing any information. Moreover the two orthogonality relations and 


:\d,mp = 0, 


(B13) 

(B14) 

(B16) 


subsequent “equal-length” relations (B4) and (B4) are the real and imaginary parts of the 


following complex orthogonality relations 


dM ■ aM = 0 dim ■ dim = o Wze 


(B16) 


where M is dehned in equation (34). They are zero on the real axis, and therefore also zero 


in the entire upper half plane. This is the key to completely solve the single-speed sector. 
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